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CALCULATION OF EQUILIBRIUM CONSTANT 
FOR DIMERIZATION OF HEAVY WATER MOLECULES 

PACS 36.40.-C, 61.20.Ja IN SATURATED VAPOR 


The magnitude and the temperature dependence of the equilibrium constant of dimerization of 
heavy water molecules in saturated vapor in terms of the second virial coefficient of the equation 
of state have been determined. An expression is found for the equilibrium dimerization constant 
of water vapor molecules, which contains terms involving the monomer-monomer, monomer- 
dimer, and dimer-dimer interaction. The obtained results are compared with experimental 
data. The equilibrium constant of dimerization in heavy water vapor is shown to exceed that 
in light water vapor within the whole temperature interval. 
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1. Introduction 

The unusual properties of water (H 2 O) have been 
known since ancient times [H-Q- a rule, they 
are explained by the existence of hydrogen bonds 
that arise between water molecules and result in the 
formation of molecular complexes, such as dimers, 
trimers, and so forth 313 . The discovery of heavy 
water (D 2 O) and its further study showed that the 
substitution of hydrogen by deuterium results in sub¬ 
stantial changes of properties in comparison with 
light water. For instance, the ternary point tempera¬ 
ture for D 2 O is by 3 K higher than the corresponding 
parameter for H 2 O, whereas its critical temperature, 
on the contrary, is by 4 K lower. The volatility of 
heavy water is lower than that of light water. Heavy 
water is more hygroscopic. Even the 30% solution of 
heavy water in light one is toxic and leads to the 
death of live organisms. When studying the proper¬ 
ties of the vapors of light and heavy waters, it is 
rather successful to apply the virial equation of state, 
which is usually confined to the second virial coef¬ 
ficient. The values of second virial coefficient found 
experimentally for light and heavy waters are dif¬ 
ferent by a factor of 2-^3 in the whole temperature 
interval Q- One should expect that this fact can 
considerably influence the value of equilibrium con- 
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stant for the molecular dimerization in heavy water 
vapor. 

Note that the differences between H 2 O and D 2 O 
manifest themselves not only in thermodynamic but 
also kinetic properties. In particular, the kinematic 
viscosity along the saturation curve of heavy water 
is about 25% higher than the kinematic viscosity of 
light water in the whole temperature interval Q. This 
work aimed at calculating the dimerization degree 
in saturated D 2 O vapor with the help of the second 
virial coefficient in the equation of state. 

2. Determination of the Dimerization 
Constant for Molecules in Saturated Vapor 

It is well known that the equilibrium properties of the 
dimerization process {m + m 4A d) are described by 
chemical thermodynamic methods. According to the 
latter, the chemical potentials of monomers, and 
dimers, p-d, satisfy the equality 


Td — ( 1 ) 

At the same time, they are functions of the corre¬ 
sponding concentrations. Therefore, Eq. (1) is actu¬ 
ally an equation for the indicated concentrations. The 
molar concentrations of water monomers and water 
molecules united in dimers are defined as 

Cm “ ^m/^0: ~ 
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respectively, where 

71q = Tlfn -j- 27l(i = ^ (3) 

is the initial density of water monomers in non- 
dimerized vapor. According to the concentration def¬ 
initions for monomers and dimers (2), the condition 
of their normalization looks like 

Cm — 1- (4) 

It was shown in work Q that, at small deviations of 
saturated vapor from the ideality, the concentration 
of dimers is determined as follows: 


Cd ~ C + 


( 5 ) 


where C = and Kp{T) is the dimeriza¬ 

tion constant. In the general case, the chemical po¬ 
tentials of components in a mixture of monomers 
and dimers contain additional contributions associ¬ 
ated with the interaction between particles, which 
look like 


(id) I (ex) 

= N +^J-i , 


( 6 ) 


where i = m,d. If the dimerization constant and, ac¬ 
cordingly, C = 2inoTKp[T) are unknown, the combi¬ 
nation of Eq. (1) with the equation of state 


parameters of the van der Waals equation of state for 
a gas mixture, which describes the excess pressure 
induced by the monomer-monomer and monomer- 
dimer attraction forces, respectively. In the quadratic 
approximation in the concentration Cd, we obtain 

C = Co + ^Co + •■•I (10) 


where 

h = P 2 + 2 pi X 

^ - (Jail - 4 qi 2 + 0-22) fr \ 

y 4 (^-l/riQ + - (ai 2 - 2aii)'j /TJ 

P2 = 2no{0.25vQ^'’ — {aii—ai2+0.25a22)/T). 

The linear approximation in the concentration Cd 
contains the contributions that involve for only the 
monomer-monomer and monomer-dimer attraction 
(the parameters an and 012 of the equation of state 
in Eq. (8)). At the same time, the quadratic approxi¬ 
mation also includes contributions which are a conse¬ 
quence of the dimer-dimer interaction (the parameter 
022 in the equation of state). For the saturated vapor 
of light water, the account for those interactions is 
crucial 0. 


P — noT(l -|- noBg.xp{T) + ...) 


(7) 3. Calculation of the Dimerization Constant 


allows one to obtain an explicit expression for 
Kp{T). In work Q, a relation between the second 
virial coefficient Be-^piT) in the equation of state, the 
dimerization constant Kp(T), and the parameters of 
the intermolecular interaction in partially dimerized 
water vapor was established. In the linear approxima¬ 
tion in the concentration Cd, we have 


C — Co, 
Co = 


Bexp(r)-u^'"^ + aii/r 


( 8 ) 


~ “ (“12 “ 2aii)/T - l/(2no)) ’ 

where 

pi = l + 2no(u^”"^ - (an - iai2)/T), (9) 

and are the excluded volumes of a monomer 
and a dimer, respectively; and an and ai 2 are the 
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In order to determine the quantity ( = 2nQTKp{T) 
or, equivalently, Kp{T), we need to know the experi¬ 
mental values of second virial coefficient Be^p{T), ex¬ 
cluded volumes (i = m,d), and gravitation con¬ 
stants amn {m,n = 1,2) in the van der Waals equa¬ 
tion. The value of second virial coefficient for the sat¬ 
urated vapor of heavy water was calculated proceed¬ 
ing from the experimental data on the pressure, den¬ 
sity, and temperature at the saturation curve . The 
quantities Vq'^ {i = m,d) and amn {m,n = 1,2) are 
connected with the behavior of intermolecular inter¬ 
action potentials. Let us take into account that wa¬ 
ter monomers and dimers are permanently rotate in 
the gaseous state, so that the microscopic potentials 
determining the interaction between water molecules 
and dimers become effectively averaged [IMl- A de¬ 
tailed discussion of the potential averaging over the 
monomer and dimer orientations can be found in 
works l4ll5l |. 
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For the averaged potentials of interaction between 
water monomers and dimers, let us take advantage of 
the Sutherland potential 


Uzj{r 12) = 



7*12 Cij , 
^12 > Cij, 


( 11 ) 


where i,j = m,d. In this case, the quantities 
and an can be found, by using the known procedure 
0- As a result, they equal 



Fig. 1. Equilibrium configuration of a water dimer 


V 


(m) 

0 


Ih^r 2 (^Yn) 

—an = EmV^ 


In all further calculations, we suppose that the ex¬ 
cluded volumes of monomers and dimers coincide 
with the four-fold volumes of hard spheres with the 
radii Vm = 1.58A and = 2.98A, respectively. In 
addition, we adopt that the averaged values of in¬ 
teraction constant equal to those quoted in Table 1 
(see work 0 ). Those values are approximately four 
times larger than the constant of dispersion interac¬ 
tion |10l4l3l |. because the dipole moment of each wa¬ 
ter molecule is induced by both the fluctuations of 
the electron density at neighbor water molecules and 
changes in the orientations of the bare dipole mo¬ 
ments of those molecules (see works 00 ). 

In order to find vf''^ and 022 , the rotation of dimers 
has to be taken into consideration. Therefore, the 
dimer radius should be taken equal to the monomer 
diameter: = a mm- The average polarizability of a 

rotating dimer ad ~ Therefore, Sd ~ Then 

(d) _ 167r ^3 o„im) _ (d) „„ (m) 

Vq - ^ ^d 5 ^22 - 


Making allowance for dimer rotations, the interac¬ 
tion between a dimer and a monomer is described by 
the parameters rdm = and £dm = 2£mm, which 
brings us to the formula 


27 (m) 

ai 2 ~ ^ ^0 ^mm- 


The results of calculations are presented in Table 2. 


4. Discussion of the Obtained 
Results and Conclusions 

The analysis of the data quoted in Table 2 and 
work @ demonstrates that the difference between the 


dimerization constants of heavy and light water va¬ 
pors is rather substantial. Depending on the temper¬ 
ature, they differ from each other by a factor of 2 
to 3. From the principal viewpoint, this is a result of 
the difference between the character of thermal exci¬ 
tations in heavy and light water dimers. Concerning 
the corresponding parameters of dimer ground states, 
they are close to one another. At the same time, the 
rotational quanta of heavy and light water dimers 
Qi = fi^/2Ji, where h is Planck’s constant, and A 


Table 1. Averaged constant of interaction 
between water molecules, €m (^b is the Boltzmann 
constant, Tc is the critical temperature) 


T 

300 K 

400 K 

500 K 

600 K 

EmlKBTc 

3.08 

3.05 

2.70 

1.78 


Table 2. Dimerization degree 

and the dimerization constant in the saturated 
vapor of heavy water 


T,K 

dd (D 2 O) 

Kp{T) (D 2 O) 

300 

0.005 

0.1580 

325 

0.016 

0.1279 

350 

0.034 

0.0858 

375 

0.064 

0.0593 

400 

0.102 

0.0387 

425 

0.149 

0.0257 

450 

0.202 

0.0166 

475 

0.251 

0.0104 

500 

0.295 

0.0063 

525 

0.343 

0.0038 

550 

0.373 

0.0022 

575 

0.439 

0.0014 

600 

0.552 

0.0009 

625 

0.886 

0.0008 
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Fig. 2. Temperature dependences of the dimerization constant 
for the saturated vapors of heavy ( 1 ) and light {!') waters 
according to the results of work IlSll and for heavy (2) and 
light { 2 ') waters according to the calculations by formula (10) 


Table 3. Rotational quanta of dimers (in cm ^ 
units) for various intermolecular interaction potentials 


Potential model 

H 2 O 

D 2 O 

Ox 

Oy 

Oz 

Ox 

Oy 

Oz 

GSD 

0.21 

8.65 

0.21 

0.18 

4.32 

0.21 

SPC 

0.23 

8.49 

0.22 

0.20 

4.24 

0.24 

SPC/E 

0.24 

8.39 

0.22 

0.20 

4.19 

0.25 

TIPS 

0.24 

9.42 

0.22 

0.20 

4.71 

0.25 

TIP3P 

0.24 

9.44 

0.22 

0.20 

4.72 

0.25 

SPCM 

0.21 

9.58 

0.23 

0.18 

4.79 

0.21 


Table 4- Frequencies (in cm~^ 

units) of vibrations in H 2 O and D 2 O dimers 

for the SPC and TIPS potentials 


SPC 

TIPS 


H 2 O 

D 2 O 

H 2 O 

D 2 O 

LJl 

70.51 

49.86 

73.71 

52.12 

UJ 2 

240.73 

170.22 

212.63 

172.24 

^3 

246.00 

219.58 

243.59 

201.72 

UJ 4 

310.53 

233.37 

322.32 

227.19 


is the moment of inertia with respect to the i-th axis, 
are different. In the case of the dimer configuration 
depicted in Fig. 1, the rotational quanta of dimers 
for various potentials of intermolecular interaction are 
indicated in Table 3. The corresponding differences 


amount to 15-20% for rotations about the axes x 
and 2 and approximately 100% for rotations about 
the axis y [l7|. 

The energies of vibrational excitations are also 
considerably different. The frequencies of small vi¬ 
brations in H2O and D2O dimers are compared in 
Table 4. 

In Fig. 2, the results of calculations of the dimeriza¬ 
tion constant for light and heavy water vapors carried 
out on the basis of the second virial coefficient in the 
equation of state are shown, as well as the results 
of direct calculations of the dimerization constant 
using the statistical physics methods by determin¬ 
ing the internal partition functions of monomers and 
dimers [13 . It is evident that the dimerization con¬ 
stants determined on the basis of experimental val¬ 
ues for the second virial coefficient correlate well with 
the results of theoretical calculations of the dimeriza¬ 
tion constants obtained in works One can see 

that the equilibrium constant of dimerization of heavy 
water molecules substantially depends on the effects 
of the interaction between monomers and dimers. It 
is owing to this interaction that the dimerization of 
molecules takes place. On the basis of the experimen¬ 
tal values of second virial coefficient, we obtained the 
value of dimerization constant. 

Attention should be paid to the fact that the tem¬ 
perature dependences of the dimerization constant for 
light and heavy waters, which were calculated us¬ 
ing different methods, have an opposite relative ar¬ 
rangement. Moreover, at temperatures in a vicinity 
of the ternary point, a considerable discrepancy is 
observed between the values obtained for Kp{T) by 
different methods. Unfortunately, we cannot explain 
now the origin of this difference, because the values 
of second virial coefficient in the temperature inter¬ 
val 300-459 K, which were used at calculations, are 
not quite reliable. It is so because I) there is a dis¬ 
crepancy in the determination of experimental values 
for the parameters in the equation of state Blii, 
and 2) there is no possibility to verify the values 
of Bexp{T) with the help of experimental values ob¬ 
tained for the viscosity of heavy water vapor. In addi¬ 
tion, in the case of the direct calculation of Kp{T), the 
determination accuracy for vibrational frequencies is 
directly connected with the choice of intermolecu¬ 
lar potentials. Unfortunately, it is difficult to specify, 
which of the potentials used in the literature is the 
most adequate. Moreover, intramolecular vibrational 
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and rotational modes are considered independent in 
work [^. However, at the large values of rotational 
quantum number, the dimer parameters considerably 
differ from their values in the ground state, which 
should be accompanied by changes in the vibrational 
and rotational modes. At temperatures higher than 
400 K, the relative influence of indicated factors de¬ 
creases, and the values of equilibrium dimerization 
constant for saturated vapors of light and heavy wa¬ 
ter molecules agree with the results of direct calcu¬ 
lations carried out in work 0 . We intend to study 
those issues elsewhere in more detail. 

The authors express their sincere gratitude to 
Prof. M.P. Malomuzh for the fruitful discussion of 
ideas and results of this work. 
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JI.A. Byjiaein, C.B. Xpanamuu, B.M. Max/iaunyK 

POSPAXYHOK KOHCTAHTH 
PIBHOBAFH ;i,HMEPI3ALJ,n MOJIEKYJI 
HACHHEHOI HAPH BAHCKOI BOJjyi 

P e 3 lo M e 

Po6oTa npHCBnuena BHSHaueHHio BejiHUHHH Ta TeMnepaxyp- 
Hoi’ 3ajie:?KHOCTi KOHCxaHTH piBHOBarn ;i;HMepH3aii;h MOJieKyji 
HacHueHOi napn BaxcKoi" bo^h Bi^noBi^no ;i,o ^pyroro Bipi- 
ajiBHoro Koe4)iii;ieHTa piBHHHHH cxany. Snaii^eHO BHpa3 
KOHCxaHTH piBHOBarH ;i;HMepH3au;h MOJieKyji Bo;],HHoi‘ napH, 
HKHH MiCTHTb ;],0;],aHKH, LLI,0 BpaXOByiOTb B3aeMO;i;ilO MOHOMep— 
MOHOMep, MOHOMep—^HMep i ;i,HMep—^HMep. HpoBe^eno nopiB- 
HHHHH oxpHMaHHX peByjibxaxiB 3 eKcnepHMenxajibHHMH ^ann- 
MH. HoKasano, uj,o y Bciii objiacxi xeMnepaxyp KOHCxaHxa piB- 
HOBarH ;i;HMepH3au;h napH BaxcKoi" bo;i;h nepeBHiqye KOHCxanxy 
piBHOBarH ;],HMepH3aLi,h MOJieKyji napn jierKoi* bo;],h. 
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